A narrow plane cut near the crystal surface increases the probability of capture into 

the stable channeling motion up to 99 percent 
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It is shown that a narrow plane cut near the crystal surface considerably increases the probability 
of capture into the stable channeling motion of positively charged particles entering a crystal at 
angles smaller than a quarter of the critical channeling angle with respect to the crystal planes. 
At smallest incidence angles the capture probability reaches 99 percent. A pair of crystals with 
cuts bent in orthogonal planes allows to reach a 99.9 percent efficiency of single-pass deflection of a 
proton beam with an ultra small divergence. Conditions necessary for efficient single-pass deflection 
of protons from the LHC beam halo are discussed. 

PACS numbers: 61.85.+p, 29.27.Ac, 41.85.Ar, 52.38.Ph, 61.80.Fe, 41.75.Ht 



The phenomenon of positively charged particle chan- 
neling in crystals has already found several applications 
in high energy physics. In particular, the deflection of 
channeling protons and ions by bent crystals is widely 
used for beam extraction from high-energy accelerators 
P, Q as well as opens up a possibility to increase the ef- 
ficiency of beam halo cleaning at the stage of LHC lumi- 
nosity upgrade Next, the radiation of channelingpar- 
ticles, in that number, in crystal undulators [II 13, H H] , 
allows to generate bright gamma-beams. Both planar 
and axial channeling can be also used for positron radia- 
tive cooling 0]. 
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ions which we will not mention further for short) most 
intensively scattered by nuclei. To change the situa- 
tion and to reach a 99% efficiency of proton capture into 
the stable channeling motion we suggest in this letter to 
use a narrow plane cut near the crystal entrance surface 
shown in Fig. 1. Similar cuts, widely used in channel- 
cut monochromators, preserve a single crystal structure 
on the cut opposite sides. In other words, atomic planes 
behind the cut are direct continuations of that ones in 
front of it. 
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FIG. I: A narrow cut between the planes z = Z\ and z = z%, 
which preserves a single crystal structure in front and behind 
itself, allows to considerably reduce the proton transverse en- 
ergy. 

The efficiency of all these applications is restricted by 
an incomplete beam capture into the regime of stable 
channeling motion, the probability of which hardly ex- 
ceeds 85% even for zero beam divergence because of the 
fast dechanneling of, at least, 15% of protons (positrons, 
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FIG. 2: Dependencies of the averaged planar potential (thin 
black curves) and transverse proton energy (thick gray lines) 
on the transverse coordinate, 1 - transverse energy reduction 
by the crystal cut region (dashed curve), 2 - dechanneling 
peak formation by the protons with ej_ — Vmax ■ Two extreme 
approximations to the potential near its maximum are shown 
on the right. 

The protons which start their transverse motion from 
the regions of high nuclear number density, such as 
|x| < 0.21 and \x - d\ < 0.2l in Figs. 2 and 3a, ef- 
ficiently dechannel at any, even zero beam divergence, 
since such protons recur to the regions of high nuclear 
density - see Fig. 3d. The cut prevents the recurrence 



by reducing transverse energies of most of the ready to 
dechannel protons. At first, passing through the crystal 
plate < z < zi protons acquire transverse velocities 
i>i = v x (zi) directed to the channel center x = d/2, 
namely v\ > at < xq < d/2 and v\ < at 
d/2 < xq < d, where xq = x(z = 0), - see Fig. 3b. Next, 
passing through the cut, they drift towards the channel 
center (see Fig. 3c). As a result, entering the crystal 
behind the cut the protons acquire smaller potential en- 
ergy V(x 2 ) than the potential energy V(xi) they lose at 
the cut entrance (here X\ 2 = x(z 12 ))- Moving further 
with the reduced transverse energy most of the protons 
are now not able to reach the regions < 0.2A and 
I a; — d\ < 0.2 A of high nuclear density (see Fig. 3e) from 
which they would dechannel in the cut absence. 

Of course, such a transverse energy reduction is possi- 
ble if the proton incidence angle $0 = «o/ u || is sufficiently 
small. Here vq = v x (z = 0) and v\\ are, respectively, 
the initial transverse and average longitudinal velocity 
components. We will retain the later in the formulas to 
make them valid both in relativistic and nonrelativistic 
cases. In order to estimate the necessary $0 value one 
can, first, conjecture that to make a transverse kinetic 
energy ev\/2 sufficiently small (say, less than a quarter of 
the planar potential amplitude V max , reserving another 
quarter for the potential energy), the angle $1 = vi/v^ 
should not exceed a half of the critical channeling an- 
gle 9 c h = \/2V max / e, where e is the total proton energy. 
Second, in order not to considerably change the angle $1, 
at least, its sign, the incidence angle i?o should be con- 
siderably (say, twice) less than $1 . Combining these two 
arguments one comes to the condition 



(1) 



of efficient dechanneling probability reduction by the cut. 

To obtain Eq. (1) more rigorously let us develop a 
simple model of the transverse energy reduction by the 
cut. Consider a proton entering a crystal with a trans- 
verse velocity vq in a point xq located in the high nuclear 
density region. Such a proton possesses a transverse en- 
ergy e± = evl/2 + V(x ) close to V max - see Fig. 2. 
After passing through the crystal plate separated by the 
cut, the proton acquires a coordinate x\ and a velocity 
V\. Entering then the cut, it loses the potential energy 
V{x\), passes freely through the cut with the constant 
velocity v\ , and acquires a transverse coordinate 



x 2 = x(z 2 ) = xi+vi ■ (t 2 - h), 
where t\ — z\/v\\ and t 2 = z 2 /v\\, and an energy 
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FIG. 3: Proton "phase space" a) at z = 0, b) at z = z\, c) at 
z = Z2 and e) at z = Z2 + ttv\\/2u) in the cut presence and d) 
at z = z 2 + 7TVu/uj in its absence. 



correctly describing the motion of most part of protons 
in the planar potential at small incidence angles. Substi- 
tuting the corresponding x-coordinate 

x\ = d/2 + clq coscj^i + vq sinuti (5) 

and velocity x-component 

Vl = —CIqLU Sm LUt 1 + VqU! COS LUt 1 (6) 

one can represent e' ± in the form of a bilinear combination 
of vq and ao with some coefficients dep ending on z\ and 
z 2 . Here a = xq — d/2 and u = \fkje are, respectively, 
the amplitude and angular frequency of the channeling 
oscillations. The coefficient by the product aoVo vanishes 
and Eq.(3) takes the form 



e 1 = e— cot cut 
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if the cut coordinates satisfy the relation 
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We find this choice of the coordinate z 2 quite appropriate 
since Eq. (7) is much simpler than the explicit form of 
Eq. (3), it satisfies the symmetry conditions e' ± (ao) = 
£j^(— do) an d e'j_(vo) = e'j_(—Vo) and reaches a minimum 
both at ao = and vq = 0. In addition, Eq. (7) allows 
to find the coordinate 



at the second entrance to the crystal at z = z 2 . 

To evaluate the cut coordinates z% t2 which minimize 
the transverse energy e' ± analytically let us use the 
parabolic approximation 
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at which the minimal value 
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of (7) is reached. Here V(ao) — ka^/2 is the initial 
potential energy and 6(ao) = (l/v\\)^2V(ao)/e is the 
proton deflection angle. Putting a = d/2 one obtains 
V(d/2) ~ V max and 8 (d/2) ~ 6> c h. Assuming that pro- 
tons quickly dechannel at e'j_ > e±dech, where s±dech is 
an effective dechanneling energy, one indeed comes to the 
condition (1) of efficient cut use at quite natural value of 

£±dech Vmax/^' 

From Eqs. (8) and (9) one easily obtains that 
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(ii) 



at 9(d/2) ~ ch and z 2 /zi ~ 4.235 at tf = O.250 ch . 
Putting fc = 8Knaa;/d 2 ~ 53.1eV/A 2 one can easily esti- 
mate that zi value increases from 100A at e = lMeV, to 
0.2^m at e = lGeV and to 17^m at the LHC energy of 
e = 7TeV. Expressing the equation of transverse motion 
in the form ev 2 d 2 x/d 2 z = —V(x)' and observing that 

zi,2 oc v\\\/e, one can easily see that the e± evolution 
from z = to z = z 2 remains the same for any energy. 
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FIG. 4: Dependencies of the proton transverse energy at z = 
(dashed) and z = 22 (solid curves) on the initial proton 
coordinate at different angles of proton incidence with respect 
to the (110) Si plane. 

The results of such an evolution simulated for three 
incidence angles $o using the Moliere approximation for 
the averaged planar potential are illustrated by Fig. 4. A 
comparison of the solid and dashed curves demonstrates 
an efficient e' ± reduction by the cut in the regions of high 
nuclear density. At ^ < c /i/4 this reduction allows 
to reach e' ± < V max /2 and, thus, to avoid fast dechan- 
neling by 98-99% of the protons while the rest 1-2% of 
them form sharp peaks at xq — x p and quickly dechannel. 
Though the origin of these " dechanneling peaks" can not 
be explained using the approximation (4), Eqs. (7) and 
(10) correctly describe the z?o dependence of the maxi- 
mum e' ± value reached outside the peaks. In particular, 
one can see that e' ± indeed reaches V max /2 ~ 12.2eV at 
i9o slightly exceeding c fc/4. 

The dechanneling peaks originate from the protons 
with Xo — x p and e± ~ £«o/2 + V(x p ) ~ V max which 
slow down their transverse motion near the potential 
maxima acquiring very low velocities v\ and coordinates 



x\ ~ nd, n — 0, ±1,±2, .., at the cur entrance. Their 
further motion with very low transverse velocities in- 
side the cut preserves the proton transverse positions 
X2 — x\ — nd near the potential maxima. As a result, 
entering the crystal behind the cut such protons acquire 
transverse energies s' ± ~ V(x 2 ) — V max corresponding to 
the dechanneling peak height. A few outermost trajecto- 
ries in Fig. 3e belong to such protons. 

A small width of the dechanneling peaks takes its ori- 
gin from the high potential curvature k' near its maxima 
due to which even a slight x deviation from the peak po- 
sitions x p leads either to the fast proton reflection from 
the potential wall at e± < V max or to the passage above 
the potential maximum at e± > V max - see Fig. 2. In 
both cases the protons quickly joint the process of trans- 
verse energy reduction illustrated by Figs. 3 b, c and 
e. 

In order to estimate the fraction of protons contribut- 
ing to the dechanneling peaks we consider an analyt- 
ical model of proton motion in the potential V(x) ~ 
V(0) — k'x 2 /2 approximating the realistic one near its 
maximum. Substituting corresponding x-coordinate 



x,\ = xq coth Ati + vq sinh Ati , 



(12) 



where A = y/k'/e, and velocity v\ = dx(t\)/dt into Eqs. 
(2) and (3) one obtains 



/ / X Tr 1 9 V, (Xo) 



(xo-x p ) 2 , (13) 



where x p = x p (v ) is the dechanneling peak coordinate 
fixed by the equation de'j_(xo) / dxo = and 



d 2 e' ± (x ) 



d 2 x 



= -i-A(t 2 -ti) 



x[X(t 2 - ii) sinh 2 Ati + 2 coth Ati sinh Ati] (14) 

is the second s' ± derivative which directly determines the 
dechanneling peak width Ax p at the level of a dechannel- 
ing energy s±dech- This width allows to evaluate the rel- 
ative number Ax p /d of protons with e±_ > e±dech which 
we will use as an estimate 
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(15) 



of the dechanneling probability. To evaluate the later one 
has to substitute the potential curvature k' the choice 
of which is somewhat ambiguous, as an insert on the 
right in Fig. 2 demonstrates. Its lower estimate k' min = 
E m ax/xi n ~ 475eV/A 2 follows from the values of the 
maximum electric field strength E max ~ 6.17GeV/cm 
and inflection point x-coordinate x,„ ~ 0.13A, while its 
maximum value k' max = d 2 V(x)/d 2 x\ x=0 ~ 846eV/A 2 is 
reached in the point x = of potential maximum. At 
£±dech = V max /2 the probability (15) is equal to 3.24% 
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at k' — k' min and to 1.24% at k' = k' max . The averaged 
value 2.24% of these estimates is in a good agreement 
with the values 2.13% and 2.59% obtained by simulations 
for the cases of, respectively, $o — and $o = 0.5/irad. 
The probability (15) will be even lower if the protons 
effectively dechannel in a relatively short crystal at higher 
e^dech values, for instance, at e_i_ de ch = 3V max /4:. Indeed, 
Eq. (15) predicts in this case that Pdech = 2.29% and 
0.91%, respectively, at k' = k' min and k' = k' max while 
the simulations give 1.13% a,t d = and 1.34% at i9 = 
0.5/irad. 

Thus, the cut reduces the transverse energy of 98-99% 
of protons of a beam with divergence (1) down to the level 
securing their stable channeling motion. The rest 1-2% of 
the protons avoid efficient energy reduction and dechan- 
nel fast. This picture manifests itself independently of 
the proton energy provided the cut coordinates which 
are, naturally, energy dependent, are chosen according 
to Eqs. (8) and (9). Thus, to illustrate the general fea- 
tures of the cut functioning by Monte Carlo simulations 
any proton energy may be chosen. Let us therefore con- 
sider the LHC proton deflection by bent crystals as the 
most challenging channeling application. 
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FIG. 5: Dependence of the 7 TeV proton dechanneling prob- 
ability in a 1 cm bent Si crystal with (solid) and without 
(dashed curve) the cut on the r.m.s. incidence angle. 

In previous deflection experiments the low capture 
probability was predetermined simultaneously by the 
high incident beam divergence, relatively thick amor- 
phous near-surface crystal layer and fast dechaneling of 
protons moving in the regions of high nuclear density. At 
low capture probability a high deflection efficiency could 
be realized only in the multi-pass mode [l|, Q. How- 
ever the divergence of the beam fraction encountering a 
crystal placed in the beam halo of a high energy proton 
accelerator can, in fact, be made very low [I,!]. The 
amorphous layer also can be made as thin as 0.01/zm [l[. 



In the circumstances the fast dechaneling becomes the 
only obstacle impeding an efficient single pass proton de- 
flection and since the cut drastically reduces the number 
of quickly dechaneling protons, its use has to make an 
efficient single pass deflection possible. 

To demonstrate this we simulated the 7TeV LHC pro- 
ton deflection by a 1 cm Si crystal with 100m bending 
radius and cut coordinates z\ — Yl^im and = 71/iTO 
evaluated using Eqs. (8) and (9). To make the sim- 
ulation process more adequate and efficient we treated 
the large-angle single Coulomb scattering separately from 
the small-angle multiple one. The angular distribution 
of protons in the incident beam was assumed to be both 
Gaussian and cylindrically symmetric. The simulated de- 
pendence of the dechanneling fraction on the root mean 
square (r.m.s.) divergence angle i9 rms is given in Fig. 5. 
On the average the cut begins to suppress the dechan- 
neling process starting from $ rTOS ~ 0.9/irad ~ O.340 c ^ 
while the highest efficiency of its use is reached at i} rms < 
0.3/irad ~ 0.19 c h when the dechanneling probability de- 
creases down to 1% as predicted by the analytical model. 




FIG. 6: A pair of 1 cm Si and Ge crystals with cuts bent in 
orthogonal planes with 100m bending radii allows to reach a 
99.9% single-pass proton deflection efficiency at the LHC. 

The simulations also show that the protons dechan- 
neled in the crystal which is bent in horizontal xz plane 
acquires a r.m.s. divergence angle of about 0.7/irad in 
the vertical yz plane. Such a divergence allows to use a 
second crystal bent in the yz plane to deflect about 90% 
of these protons, as is shown in Fig. 6. Though the high 
beam divergence impedes the dechanneling suppression 
by the cut, a combined single-pass deflection efficiency 
by a pair of Si and Ge crystals with cuts approaches 
99.9%. Both the single passage though the crystal and 
the proton moving off the high nuclear density regions 
will drastically diminish the rate of nuclear reactions ac- 
companying the proton deflection by crystals with cuts. 
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Thus, to reach a really high single-pass deflection effi- 
ciency the beam divergence angle has not to exceed much 
0.1/xrad. The angle which the proton beam axis forms 
with the crystal planes must be of the same order of 
value. Since the best goniometers allow to position a 
crystal only with a microradian accuracy, the necessary 
beam alignment along the crystal planes can most easily 
be reached by a proton deflection by a magnetic held with 
the typical integral value of only 0.01T • m which should 
be adjusted to maximize the deflection probability. 

The symmetry of the cut region from Fig. 1 prevents 
the dusturbance of e±_ reduction process by a distortion 
accompanying the cut formation. However the crystal 
bending in the horizontal plane violates the symmetry 
and can distort the cut region upsetting the £j_ reduc- 
tion process. The simplest possibility to avoid such a 
distortion is to leave the cut region straight, starting the 
crystal bending at z > z%. To prevent the distortion 
of the cut region if it is nevertheless bent we suggest to 
tightly hx a crystal plane or set of planes with random 
crystal structure orientation and elastic properties iden- 
tical to that of the removed crystal plane inside the cut. 
The method of high-resolution x-ray scattering can 
be used to check both the crystal plane matching and 
alignment on the opposite sides of the cut. 

Note briefly that the reduction of both the dechannel- 
ing probability and averaged transverse energy by the 
cut fabricated in the crystals used to generate either 
usual channeling radiation or channelled positron radia- 
tion from crystal undulators [3j, ]% |5j, la] will diminish the 
intensity of radiation from non channeled positrons, make 
the radiation spectrum more narrow, increase the effec- 
tively used crystal length and reduce the rate of positron 
interaction with nuclei. 

The positron transverse energy reduction by the cut 
also allows to increase the number of positrons involved in 
the process of radiative cooling being effective only in the 



regions of negligible nuclear density [7| ■ Note that the cut 
can be used to improve the efficiency of radiative cooling 
both in the planar and axial cases. At that the two- 
dimensional distribution of both the nuclear density and 
averaged potential allows to reach even higher reduction 
of the quickly dcchancling proton or positron fraction in 
the axial case than in the planar one. 

Note also that the cut application to the case of neg- 
atively charged particles will, on the opposite, several 
times and more increase the dechanneling and nuclear 
reaction probabilities both in planar and axial cases. 

In principle, dechanneling probability reduction by the 
cut can both increase the depth and modify the profile 
of non relativistic ion implantation. However the small 
lengths z\ 2 ~ 10 — lOOnm will seriously complicate the 
cut production. A crystalline surface layer with low av- 
erage atomic number can be used instead of such narrow 
cuts. The crystal planes of this layer also have to be the 
direct continuations of that in the crystal. A vast expe- 
rience in heterostructure growth allows one to argue 
that a high quality GaAl layer can be produced on a Ge 
crystal most readily. 

To summarize, the fabrication of a narrow plane cut 
near the crystal surface allows to ten times and more 
decrease both the dechanneling probability and nuclear 
reaction rate of positively charged particles in crystals. 
This possibility can many times increase the efficiency 
of both the beam deflection and halo cleaning at proton 
and ion colliders, of the gamma-quantum emission by 
positrons in crystal undulators and of the positron beam 
cooling. 
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